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We apply the method of matched asymptotic expansions to analyse whether cosmological vari¬ 
ations in physical ‘constants’ and scalar fields are detectable, locally, on the surface of local grav¬ 
itationally bound systems such as planets and stars, or inside virialised systems like galaxies and 
clusters. We assume spherical symmetry and derive a sufficient condition for the local time variation 
of the scalar fields that drive varying constants to track the cosmological one. We calculate a num¬ 
ber of specific examples in detail by matching the Schwarzschild spacetime to spherically symmetric 
inhomogeneous Tolman-Bondi metrics in an intermediate region by rigorously construction matched 
asymptotic expansions on cosmological and local astronomical scales which overlap in an interme¬ 
diate domain. We conclude that, independent of the details of the scalar-field theory describing 
the varying ‘constant’, the condition for cosmological variations to be measured locally is almost 
always satisfied in physically realistic situations. The proof of this statement provides a rigorous 
justification for using terrestrial experiments and solar system observations to constrain or detect 
any cosmological time variations in the traditional ‘constants’ of Nature. 
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I. INTRODUCTION 


In recent years there has been growing interest in the observational and theoretical consequences of time variations in 
the values of the traditional constants of Nature, notably of the fine structure constant, a, [1-10], the electron-proton 
mass ratio, /i = m e /m pr , [11-13], and Newton’s ‘constant’ of gravitation, G, [14]. In all cases the experimental evidence 
that can be brought to bear on the problem is a combination of local (laboratory, terrestrial, and solar system) and 
global (astronomical and cosmological) observations [15-19]. The first theoretical challenge is to develop self-consistent 
extensions of general relativity which incorporate varying ‘constants’ rigorously by including the gravitational effects 
of the variations and ensuring that energy and momentum are totally conserved by the variations that replace the 
former constants. This is achieved by regarding the ‘constant’ as a scalar field with particular couplings and self¬ 
interaction. In the case of varying G, the Brans-Dicke theory [20, 21] provides the paradigm for a scalar-tensor theory 
of this type. Recently, the same philosophy has been applied to produce simple extensions of general relativity which 
self-consistently describe the spacetime variation of a, [22, 23], and /j,, [24]. It is also possible to extend these studies 
to include simultaneous variations of several constants, and to include the weak coupling by a generalisation of the 
Weinberg-Salam theory to include spacetime varying coupling ‘constants’ [25, 26]. In these theories the analysis of 
the behaviour of their solutions is simplified because we know that the allowed variations in constants like a are 
constrained already to be ‘small’ and will not have any significant effect on the expansion dynamics of the universe 
in recent times. The Brans-Dicke theory is different. Small variations in G will always have direct consequences for 
the expansion dynamics of the universe. Typically, a power-law time variation of G oc t~ n creates a variation of the 
expansion scale factor that goes as a{t) oc t( 2-ra )/ 3 in a dust-dominated Friedmann universe [27]. 

These theories are confronted with a variety of laboratory, geochemical, and astronomical observations. In the 
case of variations of a we have laboratory constraints on atomic lines, indirect bounds from the Oklo natural reactor 
operation 1.8 billion years ago, radioactive decay products in meteoritic data back to 4.6 billion years ago, and quasar 
spectra out to redshifts z < 0(4), as the prime sources of observational evidence against which to test theories which 
permit spacetime variations [15-19]. Generally, the data from all these diverse physical scales are lumped together 
and used to test time variations of a. Thus laboratory or solar system evidence is compared directly with quasar 
data and used to constrain the allowed cosmological variations of a. Similar tactics are used to constrain the allowed 
variations of other constants, like G or pi. 

This simultaneous use of terrestrial and astronomical bounds on constants assumes implicitly the unproven require¬ 
ment that any variation of a constant on cosmological scales is ‘seen’ locally inside virialised structures like galaxies 
or solar systems, and has a measurable effect in laboratory experiments on Earth. It is not obvious a priori that this 
need be the case: we would not expect to test the expansion of the universe by measuring the expansion of the Earth. 
The central question that this paper addresses is the extent to which global variations of ‘constants’ on cosmological 
scales that take part in the Hubble expansion of the universe are seen locally on the surface of gravitationally-bound 
structures, like planets, or inside bound systems of stars like galaxies [28]. Only if we can show that cosmological 
variations have calculable local effects will it be legitimate to use laboratory and solar system observations to constrain 
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theories of varying constants in the way that is habitually done, without proof, in the literature. So far, detailed anal¬ 
yses of spatial variations of constants have been made only for small variations, where the isotropy of the microwave 
background places very strong limits on spatial variations because of the associated Sachs-Wolfe effects created by 
the gravitational potential perturbations that accompany spatial fluctuations in ‘constants’ via their associated scalar 
fields because of the coupling of the latter to matter [29, 30]. 

We are concerned with the dynamics of spacetime scalar fields that are weakly coupled to gravity and matter. 
We will not consider theories where there are two or more scalar fields interacting amongst themselves, although 
the method we use here could also be easily extended to that scenario. Theories which introduce varying constants 
self-consistently into Einstein’s conception of a gravitation theory do so by associating the ‘constant’ C to be varied 
with a new scalar field, so C —> C(y>). The variations of this scalar field gravitate and contribute to the curvature of 
spacetime like any other form of mass-energy. They must also conserve energy and momentum and so their forms are 
constrained by a covariant conservation equation for the scalar field. Typically, this results in a wave equation of the 
form 


a<p = \f(ip)L(p : p), (1) 

where <p is a scalar field associated with the variation of some ‘constant’ C via a relation C = is a dimensionless 

measure of the strength of the space-time variation of C, f(ip) is a function determined by the definition of ip, and 
L(p,p) is some linear combination of the density, p, and pressure, p, of the matter that is coupled to the field <p and 
f(ip) — 1 for small variations in and C. This form includes all the standard theories for varying constants, like G, 
a , and p, of refs. [20, 22-24], 

We shall refer to our scalar field as the ‘dilaton’, denote it by <f>, and analyse the form of the standard equation (1) 
a little further by assuming that <j)(x, t ) satisfies a conservation equation that can be decomposed into the form 

= ( 2 ) 

where T is the trace of the energy momentum tensor, T = T (with the contribution from any cosmological constant 
neglected). We absorb any dilaton-to-cosmological constant coupling into the definition of V(<j>). The dilaton to matter 
coupling B(<j)) and the self-interaction potential, V(<f>), are arbitrary functions of <f> and k = 87 iG and c = K = 1. This 
covers a wide range of theories which describe the spacetime variation of ‘constants’ of Nature, and includes Einstein- 
frame Brans-Dicke (BD) and all other, single field, scalar-tensor theories of gravity. In cosmologies that are composed 
of dust, cosmological constant and radiation it will also contain the Bekenstein-Sandvik-Barrow-Magueijo (BSBM) 
of varying a, [23], and other (single dilaton) theories which describe the variation of standard model couplings, 
[19]. Note that one could, for example, generalise the form of this conservation equation, (2), whilst maintaining 
relativistic invariance, by adding a coupling to \/T a ^T a p, or we could also break local Lorentz invariance by adding 
extra couplings to the pressure, P, defined w.r.t. to some preferred coordinate system. In this paper we will mostly be 
considering spacetimes in which the pressure of the matter vanishes, P = 0, and so all of the potential extra couplings 
mentioned above will reduce to a form that is included in the conservation equation ( 2 ) that we assume. 

We will assume that the background cosmology is isotropic and homogeneous and work with a Friedmann-Robertson- 
Walker (FRW) background metric: 


( 3 ) 

where k is the spatial curvature parameter. In the background universe the dilaton field is also assumed to be 
homogeneous (so (j> = (j> c (t)) and therefore satisfies the ordinary differential equation (ODE): 

(a 3 <j) c ^j = B t4 ,{(j> c )K (e c - 3P c ) - V# (<j > c ). (4) 

The dilaton conservation equation also reduces to a similar ODE if we are only interested in its time-independent 
mode 4>(r) in a static spherically symmetric background, like the Schwarzschild metric. In the actual universe, however, 
spacetimes that look static in some locality must still match smoothly on to the cosmological background spacetime 
on large scales. Hence, if we are to model the evolution of the dilaton field accurately in some inhomogeneous region, 
embedded in a homogeneous background, we must demand that at large distances <j> —> (/> c in some appropriate way. 
Even if spherical symmetry is assumed, the conservation equation is generically a non-linear, second-order, partial 
differential equation (PDE) and its solution is far from straightforward to find, either exactly or approximately. Even 
numerical models are technically difficult to set-up, see ref. [34], and only allow us to consider one particular choice of 
V((/>), and the spacetime background at a time. However, as we have mentioned above, if we are to bring all of 
our experimental evidence to bear on these models, we need to know how to interpret local observations in the light 


ds 2 = dt 2 -a 2 (t) 


dr 2 


1 — hr 2 


+ r 2 {d 0 2 + sin 2 6 >d 0 2 } 
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of our cosmological ones. Perhaps the most important piece of knowledge we would like to have is the correlation 
between the local and global (or cosmological) time variation of the dilaton. In particular, we would like to find the 
criteria under which it is true that 


#x,i) ~ (5) 

that is, when does the local time-variation of <j> track the cosmological one? 

In this paper we will try to answer this question by applying asymptotic methods commonly used in fluid dynamics 
in order to construct asymptotic approximations to the behaviour of <j> close to a spherical static mass, that match to 
the cosmological solution, </> c , at large distances. From this analysis, we can derive a sufficient condition for eq. (5) 
to hold. We limit ourselves in this paper to considering spacetime backgrounds that are a spherical symmetric. In a 
subsequent work we will generalise our results to deal with more general non-spherical backgrounds, [49]. 

Throughout our analysis we will refer to the spherical, static mass as a ‘star’, however it could be taken to be 
a planet (e.g. the Earth), a black-hole, or even a galaxy or cluster of galaxies. We are mostly interested in the 
(realistic) case where the surface of our ‘star’ lies far-outside its own Schwarzschild-radius. By applying our results to 
the black-hole case, however, we will comment on the problem of ‘gravitational memory’; that is, is the cosmological 
background value of </> on the horizon at the time when a black hole forms, frozen-in, or ‘remembered’, when the 
black-hole forms, or does it continue to track the background cosmological evolution? 

This paper is organised as follows: in section II we review some of the previous studies into the problem of local vs. 
global dilaton evolution and note where our work extends and improves these studies. In section III we will introduce 
the method of matched asymptotic expansions which we will use to carry out our study, and provide some simple 
examples of its application. In section IV we define our geometrical set-up of a star in a background cosmology, 
and detail our particular choices of possible spacetime backgrounds. Then, in section V, we construct overlapping 
asymptotic expansions to study the constraints required if the local and global evolution is to match together. In 
section VI we derive the conditions that must be satisfied for our method to be valid. Taking these into account, we 
interpret and generalise our results in section VII. In this section we shall also make a conjecture about a general 
condition that is sufficient for eq. (5) to hold, which will apply to more general spacetime backgrounds that those 
explicitly considered in this paper. Finally, in section VIII we use our sufficient condition to show that we do expect 
eq. (5) to hold here on Earth. 


II. PAST WORK 

Several authors have, in the past, claimed to have shown that condition (5) holds close to the surface at, r = R s , 
of a spherically-symmetric, massive body embedded, in some particular way, into an expanding universe. The results 
derived are usually only valid for a particular choice of B(tj)) and V ( <f >) and under some, usually restrictive, assumptions 
about the background distribution of matter. In this section we will review the pioneering analyses carried out by 
Wetterich in [31] and Jacobson in [35] and describe how our study will go further. 


A. Review of Wetterich’s analysis 

In ref. [31] Wetterich claimed that condition (5) would hold for any potential, V^ w \(j>), with the properties V^{(j) c ) 
B^{(j) c )Ke c . The demonstration was confined to a given background and particular potential V^ w \cj>) oc exp(— X(f>). 
It was argued that similar behaviour should be found whenever the potential dominates the cosmic evolution of the 
dilaton field i.e. B^((f> c )Ke c . In what follows we will show that condition (5) also holds in many situations 

where the potential does not dominate the cosmic evolution of the dilaton field. It is our belief, however, that the 
reasoning given in ref. [31] is incomplete and does not show that condition (5) holds even under the restrictive 
conditions specified there. We shall briefly outline the arguments given in ref, [31] below, and show where we believe 
they fail. 

Wetterich considered a universe filled with pressureless dust and a cosmological constant. Under these assumptions 
we have T = edust , the dust density, and: 


= B^(<j>)Ke dust - V y (<f>) , 

where e is the local matter density and it was assumed that the potential is 

V {w) (0) = ue~+, 
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TABLE I: Density distribution in Wetterich’s model 


Region 

Range 

€-dust 

Description 

a 

r < R s 

e = €e 

local planet 

b 

R s < r < r c 

e = 0 

intermediate space 

c 

r > r c 

e = e c {t) 

Hubble flow 


where u> ~ O • The function B ^(0) represents the coupling of the dilaton field to the dust. Experimental bounds 

on the largest allowed violations of the Weak Equivalence Principle (WEP) that will be created by a dilaton held that 
couples to the electromagnetic energy of matter suggest that |B )fl j(0)| < IIP 4 . Wetterich considered a particularly 
simple example of a spherically symmetric massive body superimposed onto the cosmological background; table I 
shows the local energy density budget in this model. For r r c , 0 takes it cosmological value 0 c (f), which satisfies: 

4>c + 3H(j) c = we“^ + B^{4> c )Ke c . 

Wetterich only considered the case where u>e~^ c B i<? j(0 c )fte c . Assuming that the scale factor a oc t n , we find: 

</>c(t) = (/> o + 21n (t/t 0 ), 

where t = to today, and 0(fo) = 0o = — ln(6n + 2)/uto) « 140 + In is the present value of the scalar held. 

For r ~ R s , dilaton held in ref. [31] is written as: 

0 = + 4>i(r,t) + (j> e {r,t), 

where this dehnes 0 e and <f>i(r,t) is the ‘local’, quasi-static held conhguration satisfying: 

V 2 0j - /. i 2 (j)i = -B^ (0 C ) ne E 0 (R s - r), 

with V 2 the 3-D Laplacian and 9 ( R s — r) the Heaviside function; /z 2 = V^(<j> c ) and /. i 2 R 2 <C 1. Near r = R s it was 

found that (j)i/<j>c ~ B i<? j(0 c )2M/B s <C 1. If 0 e (r, f)/0 c (f) <C 1 for r ~ R s then the dilaton held will satisfy condition 
(5). If \\3<j> e /\3<j> c \ <C 1 locally then, as stated in ref. [31] we will have the required result. Wetterich argues that this 
is indeed the case in his paper. However, this does not appear to be the case. Assuming 2 m/R s <C 1, from r = R s 
spacetime is approximately Minkowski and so: 

□0 e = <j> e - V 2 0 e « -/z 2 0 e - 02 + 3H0> C - J3 0(^c)e c + B^^^esOiRs - r) {(fn + <j> e ) , (6) 

□0c = tc = -3Hj> c -V$’\<f> c ). (7) 

Now, 0; <C 0c and A* 2 A 2 <C 1 and so, near the surface of our body, the /z 2 0 e term represents a negligible correction to 
the 0 e dynamics. Therefore, if D0 C D0 e we must have: 

- B i? j(0 c )e c + B t ^{(j> c )tE (<t>i + 0e)) < 0c- 

However, under Wetterich’s assumptions that B,<^(0 C ) ~ O (B^(0 C )) and e E ~ 10 3 °e c , as is the case for the density 
of the Earth, we hnd: 

B,00(0 c )e£:0; ^ 22 f MM] ^ 171 ~ IQ 13 ( ^.<i > (0 c ) \ ^ 1 

0c ~ V 10- 4 ) Rs ~ l 10~ 4 ) 

where we have taken ~ 10 -9 in the hnal deduction. Hence, we have shown that in general the condition 
[□0e/n0 c | 1 does not hold; indeed D0 e D0 C . This result is opposite to the one found in ref. [31]. We 

conclude then that Wetterich’s analysis does not prove that condition (5) holds. If we approximate our local solution 
by 0 = (j)i +0c then we have seen that correction terms, 0 e , to this solution vary on scales much smaller than 1/0 C . As 

a result we cannot conclude that 0 ss 0 C . In fact, at the epoch t = to, the asymptotic expansion of the local solution 
should be correctly written as: 

0 ~ 0o + 02( r > to), 

We shall see later by more detailed methods that, even though the above analysis fails to show it, that we should 
expect condition (5) to hold for this set-up. 
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B. Jacobson’s result 


In ref. [35], the problem of gravitational memory [32] is considered in the context of Brans-Dicke (varying-G) theory. 
If Newton’s constant, G, can and does vary over time and space then one must ask which value of G is appropriate 
for use on the horizon of a black hole after it forms. One motivation was to discover if the black hole possesses a 
type of gravitational memory, freezing-in the value of G that existed cosmologically at the moment when it formed 
in the early universe, or whether the value of Newton’s ‘constant’ on the horizon changes over time so as to track 
its changing cosmological value in the background universe. This has been investigated by several different methods 
and found to be a small effect [33] but Jacobson’s argument was that if the cosmological variation in G, and the 
related Brans-Dicke field, 0 oc G” 1 , is slow over time-scales of the order of the intrinsic length scale of the black-hole 
(~ 2 GM) then, at each epoch t = to, one can expand the cosmological value of 0 as a Taylor series in (T = t — to) 
and, to a good approximation drop all 0(T 2 ) and higher-order terms, so 

4>c ~ 0i = 0c (^ o) + <j> c (to)T. 

Jacobson noted that (j>\ is a complementary solution to the Brans-Dicke field conservation equation in an (empty) 
Schwarzschild background: 

□«0 i = 0 , 


where D s is the d’Alembertian operator for the Schwarzschild metric. Therefore (j>\ can be added to any known 
Schwarzscliild-background solution of Brans-Dicke field equations to gain a new solution. Jacobson took T to be the 
‘curvature’ (or Schwarzschild) time so that the metric is given by: 

ds 2 = (1 - 


2m 


R 


2m 


R 


-l 


— dT 2 - 1-d R 2 - R 2 {dd 2 + sin 2 9(\(jf}. 


This time coordinate diverges as we move towards the horizon, hence also diverges in this limit. The unique 

static solution of the dilaton conservation equation, which vanishes as r —* oo, is given by: 


0 2 = In 



Whilst both cf >i and 02 diverge on the horizon, Jacobson found that there is a unique linear combination of them that 
is well-defined as r —> 2m: 


4 >jac = 0 i(£) + 2 rn 0 c (fo) 02(f) = <frc(v - r - 2 mln(r/ 2 m)) ( 8 ) 

where v = t + r + 2m\n(r/2m — 1) is the advanced time coordinate. By construction, it can then be shown that 
there exists a unique solution for 0 in the Schwarzschild background that is non-singular on the horizon and has 
0j ac (r = oo, t) « 0 C (£) for the particular case where 0 c (f) oc £ oc G _1 . If the approximation used is valid, then this 
suggests that the value of G on the horizon lags slightly behind the cosmological value, but that over time-scales much 
larger than rn there is no gravitational memory when G oc £ _1 falls in this extreme fashion. 

Jacobson’s result assumes that the cosmological region can be taken to be infinitly far away, so that the entire 
spacetime is Schwarzschild; in reality the cosmological matter will become gravitational dominant over the black-hole 
at some finite-distance. Another issue with Jacobson’s method is that, since T diverges on the horizon, the expansion 
of 0 C as a Taylor series in small T will not be valid near the horizon. In realistic models the space surrounding the 
black hole will also not be totally empty, indeed there will be an accretion disk surrounding the black-hole. The time- 
scale for matter to fall into the black-hole is ordery / r 3 /2m and this is relatively short compared with the time-scale 
over which cosmological expansion occurs. Accretion of matter into the black-hole, therefore, might well result in a 
significant difference in the time-variation of G close to the horizon or cosmologically. In this paper we will consider a 
more realistic embedding of a Schwarzschild mass in an expanding universe; and build our asymptotic expansions in 
such a way that they are well-defined on the horizon. We will conclude that Jacobson’s result [35] (see also refs. [33, 34] 
for similar conclusions arrived at by different methods) does indeed give the correction behaviour of 0 whenever the 
energy-density in the region surrounding the black-hole is low enough and for more general time-variations of G(f) in 
the background universe. 


III. MATCHED ASYMPTOTIC EXPANSIONS 

The method of matched asymptotic expansions that we will use in this paper was first developed to solve systems of 
PDEs that involve multiple length scales. It is often used in the field of fluid dynamics to study systems where there is 
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thin boundary layer, inside which the length scale of variation is much smaller than outside it. It also has applications 
in the study of slender bodies, with widths much smaller than their lengths. Other problems with multiple length 
scales include the interaction of greatly separated particles and the evaluation of the influence of a slowly changing 
background field on the dynamics of a small body. The books by Cole, [36], and Hinch, [37], provide a more detailed 
treatment of this subject. In this section we will briefly introduce the method and give some simple examples of its 
applications in order to fix ideas. 


A. Asymptotic expansions 


Critical to this method is the requirement that the approximations we will work with are asymptotic expansions 
rather than convergent (Taylor-like) power-series. In general, an asymptotic expansion will not be convergent. Thus, 
we define asymptotic approximations and expansions as follows: 

Definition III.l. ]T]' Aj f fn(S) is an asymptotic approximation to f(5) as 5 —> 0, if for each M < N the remainder 
term is much smaller than the last included term: 


f(S)~E M fr 


( 6 ) 


f M (S) 


0 as 6 —> 0. 


One then writes: 


N 

m ~ 5] m as <5 —► 0. 

Definition III.2. If definition III.l holds, in principle, for all N, i.e. we can take N = oo, then the we deem the 
approximation to be an asymptotic expansion of f(8): 

m^TTMS) as <5 —> 0. (9) 

In many cases f n (8) oc (S p ) n for some constant p and we will have an asymptotic power series. It is also quite 
common to find fi(5) oc 6 J ln(d), for some i and j. The sum in eq. (9) is a formal one, and it is not required that it 
does converge. The property required by definition III.l is, however, more useful that convergence in many cases; it 
ensures that one only needs the first few terms of the expansion to create a good numerical approximation to f(5). 


B. Singular problems and ones with multiple scales 

Consider a differential operator £ x (5) which defines some function f(x,S) by C x (5)f(x,5) = 0. We can attempt to 
solve this equation by making an asymptotic expansion of f(x, S) and solving the resultant equation order-by-order 
in 5: 

OO 

f(x,S) ~ ^ f n (x)'y n (5)as 5 -> 0, 

with x fixed. A ’singular’ problem is one where the above asymptotic expansion is not uniformly valid , i.e. it breaks 
down for certain ranges of x, typically x = O (5) or x = 0(1/5). Singular behaviour such as this can be divided 
into two distinct classes. In the first, the highest derivative in C x is multiplied by some power of 5 and so can be 
ignored everywhere apart from in those regions where the variation in f(x, 5) is fast enough to ensure that the highest 
derivative make a significant contribution. In the second class, the problem has more than one intrinsic length (or 
time) scale, one much larger than the other. The application of this method to physical problems generally tends to 
fall into this latter class. In both cases one proceeds by constructing two (or more) asymptotic approximations to the 
solutions which are valid for different ranges of x, e.g. for x ~ 0(1) and x/6 = £ ~ 0(1), with 

f(x, 5) ~ fn(x)S n as 5 -> 0, X fixed, (10) 

f(x,6) ~ Y,n=o9n(08n as 5 -» 0, £ = x/6 fixed, (11) 

and solving £ x (5)f(x,5) order by order in 6 for both expansions w.r.t. to some boundary conditions. We will call 
expansion (10) the outer solution, and (11) the inner solution. The inner expansion is not uniformly valid in the 
region £ = 0(1/5), as the outer one is not valid where x = 0(5). Because of these restrictions on the size of x, we 
will only be able to apply a subset of the boundary conditions to each expansion; in general we will, therefore, be left 
with unknown coefficients in our asymptotic approximations. This ambiguity can be lifted by matching the inner and 
outer solutions in an intermediate region where they are both valid. 
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C. Matching 


We can match the inner and outer solutions together if there exists some intermediate range of x where both (11) 
and (10) are valid asymptotic approximations to f(x,5). By the uniqueness properties of asymptotic expansions, see 
[37] for a proof, if they are both valid in some region then the two expansions must be equal. We take x to scale as 
some intermediate function, 77 (c)), with magnitude between 1 and 6, e.g. r/(S ) = S a where 0 < a < 1 , and define a new 
X coordinate appropriate for the intermediate region by 


V - */-)(<) = (^y) f- 


We then write both the inner and outer approximations in terms of X and, keeping X fixed, make an asymptotic 
expansion of each of them in the limit 5 —> 0; this is called the intermediate limit. We must be careful to neglect any 
terms in the approximations, that we find in this way, which would be smaller than the first excluded term; equality 
is then demanded between the remaining terms: 


Etc fn(x(X, v ))6 n ~ Y,LA ut \x) ln {5) 

ELo 9n{Z{X,ri))6 n ~ ZL 0 hn n) (X) ln (6) 

h ( ° ut) (x) = hi: n \x) 


as <5 —^ 0; X fixed, 
as <5 —7 0; X fixed, 
for 0 < n < S. 


( 12 ) 


Usually, the precise form we choose for 77 (c)) is unimportant and the matching can be done regardless; there are some 
cases, however, where we must be more careful (see Hinch [37] for more). 


D. Simple Examples 

1. Example 1: Simple matching 


We will start with a boundary layer example where the exact solution is known. Consider y{x, 6) which satisfies: 

2Sy"+(l + 26)y' + y = 0, 

in 0 < x < 1, with boundary conditions y = 0 at x = 0 and y = e -1 at x = 1 and where ' = d/da;. The exact solution 
is: 


y{x,S) = e 



_ e ~x/2S 
_ g —1/2(5 


(13) 


If we apply the method of matched asymptotic expansions to this, we would find that the outer approximation is 
given (to all orders) by: 


y(x, 6) ~ e x ^ 0 fixed. 

We have been able to apply the boundary condition at x = 1; the x = 0 boundary condition cannot be reached given 
our assumptions about the size of x. For the inner approximation we keep £ = x/5 fixed. In terms of £ our differential 
equation reads: 


2 y,# + (1 + 2S)y^ + Sy = 0. 


The inner solution is found to be: 

y(x, S) ~ a 0 (l - e _4/2 ) + S (a,! (l - e~ s/2 ^j - a 0 c) + S 2 (a 2 (l - e~ i/2 ^j - a + O (<5 3 ), 

with £ = x/5 held fixed. Again we can only apply one boundary condition, this time the one at x = 0. The boundary 
condition at x = 1 is beyond our reach in this approximation. In this problem the choice of the intermediate scale 
turns out to be unimportant; we take 77 = 5 1 / 2 and define X = S~ 1 ^ 2 x = d 1 / 2 ^. The intermediate limit of the inner 
approximation is: 


inner approx. ~ a$ — S 1/2 a 0 X + S( ai + \a 0 X 2 ) + O (<5 3/2 ) 
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FIG. 1: The leading order, inner and outer approximations to the solution of problem given in section HID 1, with 5 = 0.00001. 
The exact solution is not visible on this plot since it lies underneath one or other of the approximations everywhere. 


The first excluded term is order <5 3 / 2 ; we have neglected all terms of, or below, this order including the exponentially 
small e ~ x / 2<5 terms. Expanding the outer approximation we have: 

outer approx. ~ 1 - 5 1 > 2 X + 5\X 2 + O (d 3 / 2 ) . 

By our matching criteria, eq. 12, we must have: 


a 0 = 1, 

cti = 0. 

If we were to perform this process to all orders we would find that all a,; = 0 for i > 1. The fully specified inner 
approximation is therefore 


n —0 

This is precisely what we would have found by performing a Taylor series expansion of the exact solution, (13), in 
the inner limit and then dropped all exponentially small terms i.e. By performing the matching we have been 

able to lift the ambiguity in the coefficients, a t , and fully specify the inner approximation. 


2. Example 2: Scalar field in an expanding universe 

Consider a spherically-symmetric, scalar field, <f>(r,t), like the dilaton, in a flat FRW cosmology with metric 

ds 2 = dt 2 - a 2 (t) (dr 2 + r 2 (dtf 2 + sin 2 9d</ 2 )) . 

To make contact with our problem we will assume that the only source term in the /> evolution equation is homogeneous 
and proportional to the background dust density, so 

(a 3 m- (rW A H 2 

9 a 3 (t) a 2 (t)r a 3 (t ) 

where ' = d/dt and ' = d/dr and with Hq the Hubble parameter at some arbitrary time t = to and A is a constant 
parameter of the theory. We simplify the problem further by considering only the matter era, where a(t) oc t 2 ^ 3 . 
As boundary conditions we take 4> —> (j>o(t) as r —> oo and </>'/a{t )| r=ro / a ( t ) = p2m/r\ = const ; /r is a dimensionless 
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constant, and m is a length scale such that 6 = 2mH 0 <C 1. We assume that ro/a(t 0 ) ~ 0(2m). We could visualise 
this problem arising from the presence of a massive body in the region a{t)r < r 0 that creates a spatial gradient in <j>. 
Near a(t)r = ro we have 0' <f >. We take our problem to be similar to one we will consider later and determine the 
relation between 0 o and 0(ro/a(t),t). 

The small parameter in this problem is 8 = 2mHo- We work at the epoch when t = to; for simplicity we take 
a[t o) = 1. In the inner approximation we define 2= a(t)r , and £o := ro/2m ~ 0(1). We assume that the 
inner solution is quasi-static, so it depends on time only through the slowly increasing cosmological t, rather than 
(t — to)/2m. We define a dimensionless time coordinate by r = Hot, r 0 = H 0 to = 2/3. In this interior region 
0 = 0i(£, t; 8) which satisfies 


±81^) = 6 2 


^2 dr ( T 2 8 t 0I ) 


^ + ^ 2 ^^ + ^2 (? 9 «) 2 ^ 


Solving this to order 0 (<5 2 ) we find: 


<Mt£) 


0 o ( r ) 


8 

S, 


S 2 




9 (t) (e- 2 e 0 ) \ 
^ ) 


+ 0 (<5 4 ) 


where 0 q (t) and (r) are ‘constants’ of integration, to be determined via the matching procedure. g(r) = 
-pzd T ( T 2 d T (f>o ) — In the exterior we define p = Hqv to be our dimensionless radial coordinate, and 0 = 0e■ 
Our boundary conditions determine the leading-order exterior term to be <J>q (t) (= 0o(t) with some abuse of nota¬ 
tion). The sub-leading order terms, 8 l 0$ are then given by O0^) = 0. So, to order 8, 


0e{t,p) ~ (j) 0 (r) + <5 


fZd7T 7 (T)X 7 (p) 

P 


0 { 8 2 ) 


where 


X 7 (p) = A{(5) cos {j3p) + B(/3) sin ((3p) where 7 = — f3 2 < 0, (14) 

X 7 (p ) = A(0) where 7 = 0, 

X 7 (p) = C(a)e~ ap where 7 = a 2 > 0, 

t 2 T 7 (t),tt + 2tT 7 (t) iT = 7 t 2 / 3 T 7 (t), (15) 


and A(/3), B(/3), C(a) are all to be determined by the matching procedure. As in the previous example, the precise 
position of the intermediate region is not important. We choose an intermediate coordinate z = 8~ 1 / 2 p = 6 1 / 2 £ and 
take the intermediate limit of both the interior and exterior approximations. By equating our two approximations in 
the intermediate region we find B((3) = C{(3) = 0 from the 0 ( 8 ) terms and 


0 ( 1 ): 

0oi T ) = 0o(r) -> g(r) = 0 , 

(16) 

O (A 1 / 2 ) : 

J™ df3A(J3)T_fp (r) = -03 

(17) 

O (S 3 / 2 ) : 

J™ d f3[3 2 A(P)T_r (r) = ^3 . 

(18) 


By differentiating twice and applying (14) we can check that conditions (17) and (18) can be simultaneous satisfied 
for some choice of A(/3). The T 7 (r) can be made orthonormal w.r.t. to some inner product and so eq. (17) can, in 
principle, be inverted to find A{/3). We omit this step, however, since we are mostly concerned with the effect of the 
exterior on the behaviour of 0 in the interior rather than vice versa. By finding the interior solution to 0 (<5 4 ) and 
the exterior to 0 ( 8 2 ) we can show that 0[ (r) = 0 .; we have now fully specified the interior solution to 0 (S 2 ): 


We have that: 


0 / (t, R = a(t)r) ~ 0o (t) - ^ ^ (r + + 0 (<5 4 ) 


0l\ar=r o _ ^ _ 16TO/Zr 0 


00 


9t 3 0O 
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and (f> o ~ 4A/9f for large t. Therefore we have shown that for the time variation of <j> at ar = ro to track the 
cosmological time variation we need 


It is clear that whatever the value of A, the rate of time variation in <f> will tend to homogeneity as t —> oo. This 
example is a greatly simplified version of the problem we will consider in the rest of this paper. 


E. Application to General Relativity 

The use of matched asymptotic expansions in general relativity was pioneered by Burke and Thorne [38], Burke 
[39], and D’Eath [40, 41] in the 1970s. These authors used them to the study how the laws of motion of a test body 
were affected by the background spacetime. We shall now outline how matched asymptotic expansions are applied in 
general relativity. 

We assume that the universe 


C=(M,g ab ,T ab ,^), 

with T ab is the energy-momentum tensor and <f> the dilaton, can be viewed as a background cosmology, 

Co= {M 0 ,g° ab ,TZ b ,<l> 0 ), 

onto which some localised, interior configuration has been superimposed in a non-linear fashion. In what follows, for 
simplicity, we will require the interior configuration to be static in some coordinate system. We demand that the 
‘size’ of the interior region be given by a single parameter m; and that as m —> 0, with { x fixed, the interior region 
disappears and C —> Cq. We shall formalise this statement later. 


1. Length scales 

The length scale of the interior region, as defined by its curvature invariants is denoted by Lj (m), with Lj (m ) —> 0 
as m —> 0. The length scale of the background (exterior) region is written Le■ For the asymptotic expansion method 
to be viable we need 5 = Lj/L# <C 1. The effect of the interior on the background configuration can then be treated 
as a linear perturbation to Co, with 5 playing the role of a small parameter. We can similarly treat the effect of the 
background on the interior as a linear perturbation. 


2. Five-dimensional manifold 

For each m, in some interval [0, m max ), we write the global configuration as C m = (M m , g ab (m), T ab (m), <j>(m )). 
Following Geroch [42] and D’Eath [40], we consider a five-dimensional manifold with boundary, A f, that is built up 
from spacetimes (M m , g ab {m)) for m £ [0,m max ); g ab { 0) = . As D’Eath noted we should properly exclude from 

M 0 a smooth time-like world line, Zq, that corresponds to the ‘position’ of the interior region; this is illustrated in 
figure 2. We require that the contravariant metrics g ab (m) define a smooth tensor held on A/"; in addition, we require 
that the dilaton, < p(m ), defines a smooth scalar held on A f. These conditions are required by our assumption that the 
interior region has only a small perturbing effect over length scales >> Lj(m). 

In an open subset on A7, we choose a coordinate chart (f, r, 9 , </>, m) such that the aforementioned world line lo is 
given by (r = 0, m = 0). In the limit (tL^ , rLf^ 1 , 6, <f> ) — > constants and m — > 0 we can give gab(tn), T ab (m) and <p(m) 
as asymptotic expansions about , T^ b and , respectively. This is the exterior approximation, and is appropriate 
for considering the perturbing effect that the interior has on the background universe. 

For each epoch t = to, we shall also dehne an interior approximation. For this, we take the asymptotic expansions of 
g ab {jn), T ab (m) and </>(m) in the limit ((t — t 0 )LJ 1 ( m), rLj 1 (m), 9, </>) —> constants, and m —> 0. This approximation 
is appropriate for considering the perturbations produced in the interior region by the background cosmology. This 
is the problem that we are most interested in. We shall assume that we know the leading-order configuration in this 
limit and in what follows we take it to be a Schwarzschild metric with a quasi-static dilaton held. 
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N 


exterior region interior region 



FIG. 2: The 5-D manifold J\f is built up from spacetimes (Mm, g a b(m)) for m £ [0 ,m ma .x)- 


3. Intermediate region and matching 

We match the interior and exterior approximations in some intermediate region, where ((f — to )Ll nt , rL~^ t ) —► const 
as m —> 0, and where Li nt = LeS 1 ^°‘ = Lj5~ a , with 0 < a < 1. In many cases the precise value of a is not important 
(see section IIIC). Following D’Eath, [40], we assume that all the functions in our approximations are sufficiently 
well-behaved in this region so as to admit a power-series expansion in the radial coordinate; given this we shall not 
need to examine this region explicitly. A typical term in the interior expansion will, in this matched region, look like 
fij(to; L~^ t (t — to)) where £ = rLJ 1 and, i and j are rational numbers; % > 0. In the exterior region, a typical 
term will look like S k p l gki(to', (t — to)), again with k and l rational and k > 0 and p = L~j^r. Matching will require 
that g(fc+i)j = fij■ When we consider the Tolman-Bondi class of background metrics we will modify this procedure 
slightly. For these metrics there are two choices for radial coordinate: 

• the physical radial coordinate R defined by surfaces of ( t , R) = const with surface area 47 tI? 2 , 

• the radial coordinate which is constant on the world-lines of dust particles, r. There is an ambiguity in this 
definition of r in that the defining property is also satisfied by any arbitrary function h(r) of r. We lift this 
ambiguity by demanding that R = r at t = to- At later times the Einstein equations give us R = R(r, t). 

The functions in our interior expansion will generally be quasi-static when written as functions of R and r. We will 
choose to expand the interior and exterior approximations as power series in the physical radius, R, when taking the 
intermediate limit. 


IV. GEOMETRICAL SET-UP 
A. General Picture 

We shall assume that the dilaton held is only weakly coupled to gravity, and so its energy density has a negligible 
effect on the background spacetime geometry. If we ignore this back-reaction then we can, rather than solving the 
full set of Einstein-matter-dilaton equations, simply consider the dilaton Held’s evolution on a fixed background with 
given matter density. We will take the background to be some known exact solution to Einstein’s equations with 
matter possessing the following properties: 

• The metric is approximately Schwarzschild, with mass m, inside some closed region of spacetime outside a 
surface at r = R a . The metric for r < R s is left unspecified. 

• Asymptotically, the metric must approach FRW and the whole spacetime should tend to pure FRW in the limit 
to — > 0 . 
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• As the local inhomogeneous energy density of asymptotically FRW spacetimes tends to zero, the spacetime 
metric exterior to r = R s must tend to a Schwarzschilcl metric with mass m . 


cosmological region 




FIG. 3: We are concerned with the evolution of the dilaton field, cj >, in spacetimes where a mass m, of radius R 3 , has been 
embedded into the cosmological fluid. 

This list of requirements is far from exhaustive, and they shall be re-expressed in more rigorous fashion later. In 
the local, approximately Schwarzschilcl, region the intrinsic length scale, Lj, of a sphere centred on the Schwarzschilcl 
mass with surface area 4 irR is given by considering the curvature invariant: 

3/2 

Lt = (±R abcd R abcd y 1/4 = (19) 

(2m) ' 

In the asymptotically FRW region, the intrinsic length scale is proportional to the inverse root of the local energy 
density: l/V ne + A. We shall assume that, in line with current observations, that the FRW region is approximately 
flat, and so we set our exterior length scale appropriate for the epoch at t = to is the inverse Hubble parameter at 
that time: 


L e = 1/H 0 . 


Our small parameter is then defined to be 


S = Li/Le- 

Formally we require that, as S —> 0, our choice of spacetime background should be FRW at zeroth-order in the 
exterior limit and Schwarzschilcl to lowest order in the interior limit. In section VI we will say more about what is 
required of the background for this matching procedure to be valid. In this paper we shall, in addition to the criteria 
given above, restrict ourselves to the subcase where the spacetime is spherically symmetric. All spherically-symmetric 
solutions to Einstein’s equations with matter, where the matter is pressureless dust and a cosmological constant and 
the flow lines of the matter particles are geodesic, fall into the Tolman-Bondi class of metrics (for a review of these 
and other inhomogeneous spherically symmetric metrics see ref [43]). This class of possible Tolman-Bondi spacetimes 
is parametrised by two arbitrary functions of one spatial variable, r. We will only look at the flat, Tolman-Bondi 
models with non-simultaneous initial singularities (which where rediscovered by Gautreau in [43, 44]), and the non-flat 
Tolman-Bondi models with simultaneous initial singularities. These two subcases are fully specified by prescribing 
only a single function: the matter density on some initial space-like hypersurface. Another metric that has often been 
used to study the effect of the universe’s expansion on solar system dynamics is the McVittie metric [43, 45]. In the 
exterior limit the McVittie metric asymptotes to a dust + A FRW cosmology. For radii r where r ^ 2m, r <C H^ 1 , 
the McVittie metric looks like Schwarzschilcl spacetime; the horizon in the McVittie metric possesses a curvature 
singularity however, and it cannot therefore be used to model a black hole in an expanding spacetime. 
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B. Case I: The McVittie background 


The earliest studies of the gravitational field produced by a spherically-symmetric, massive body in an expanding 
universe were based upon McVittie’s solution to Einstein’s equations [43, 45]. McVittie’s solution is a superposition 
of the Schwarzschild metric and a FRW background. In isotropic coordinates the metric, is 


ds 2 


I ~ n(t, r ) 

1 + M*> r ) 


dt 2 


[1 + fi(t, r)] 4 
(1 + \kr 2 ) 2 


a 2 (t) [dr 2 + r 2 {d6> 2 + 


sin 2 0d(j> 2 }\ , 


( 20 ) 


where 

Mt ’ r) = 2 ^ (l + ^ r2)1/2 ’ (21) 

and m and k are arbitrary constants. In the limit a = 1, m is the Schwarzschild mass, k is the curvature of the 
surfaces (t,r) = const in the m = 0 (i.e. FRW) limit. As a model of a black hole in an expanding universe it has 
the distinct disadvantage that the ‘horizon’ is a naked curvature singularity. This defect aside, the McVittie metric is 
believed to be a good approximation to the exterior metric of a massive spherical body with a physical radius much 
larger than its Schwarzschild radius. 

In the flat (k = 0) case, the local energy density depends only on time and is the same as the FRW energy density 
to which the solution matches smoothly as r —> oo. The pressure, however, is not the same as the FRW pressure. 
Apart from the vacuum-energy dominated case (i.e. where e = — P), it is not possible to prescribe a dust, P = 0, or 
any other, non-vacuum, barotropic, P = P(e), equation of state to hold apart from in the asymptotic FRW limit. 


1. Exterior Metric 


We shall work in the (r, t) coordinates introduced above. In the exterior dr ~ H 0 1 ~ dt. We define dimensionless 
exterior coordinates r = H 0 t and p = H 0 - Then, with 5 = Li/Le = H 0 R^ 2 (2m) 1 ^ 2 <lwe have: 


ds 2 

ext 


Hk 


{g^+Sg^ + 0(d 2 )}dx a dx b 


= H n 


1 - p) 

1 + Sp,(r,p)\ 


2 [1 + Sp(T,p)] 4 2/ X rj„2 , 2rjfl2 , 0 :„2/,n,2l 


dt — 


(1+ 1W 


( 22 ) 


■a 2 (r) [dp 2 + p 2 {d0 2 + sin 2 0d(j) 2 }] > , 


where p(r, p) = j (2 m/R s ) 3 ^ 2 {pd(r)) 1 (1 + if IkP 2 ), a(r) := a(t) and 0° = kH 0 2 . 


2. Interior Metric 


In the interior it is most convenient to work with 


R=[l + p] (1 + \kr 2 )- 1 > 2 a{t)r 

as the radial coordinate. If k = 0, then the surface (t, R) = const has area 47ri? 2 . When t = to the interior 
geometry will be approximately Schwarzschild over scales where dl? ~ 0(R S ) and d(t — to) ~ 0(Li). We define new 
dimensionless coordinates by 


T = (t — t 0 ) /Li and £ = R/R s . 

In these coordinates the interior metric is: 

d*L ~ Rl {j ( ab (0 + dj ( ab (ST, O+O (S 2 )} dx a dx b 

= ^(frM^+fW + s in 2 M0 2 }]|, 


(23) 
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where A(£) = 1 — 2m/(R s £), 

/ e» \ 1 / 2 

r/> = d£-Sh(5T)t^J A(0 1/2 dT, 

and h(ST) = a'(r = 8T)/a(r), flk(ST) = 12° a _2 (r). We have defined X(£) = | (l + A(() 1 / 2 ') 2 . Note that the interior 
metric functions j® b , j^ b etc. can, at each order, be written in quasi-static form, as functions only of £ and ST. 


3. Energy-Momentum Tensor 


The dilaton field theories we are considering couple to the matter sector through the trace of the energy momentum 
tensor, nTff = k(c — 3 P). In the exterior expansion: 


ij 0 - 2 «r; = 4(T) + 6ei(T, P ) + o(6 2 ) 


= 12/i 2 (t) + 6/i'(t) 


In the interior expansion we have: 

I? 2 «t; = S 2 e\{ST,0 


1 + Sjl(T , p) 
1 - p) 


3fife (r) 


2 - £/2(t, p) 
1 - (5/1 (r, p) 


(l + <5/z(-r,p)) ' -4^ 


= 128 


2 m 

i?7 


h 2 (ST) + 6 S 2 h'(8T)A(Z)-V 2 


When rewritten in quasi-static form, the only term in the interior expansion of T is at 0(S 2 ). 


(24) 


(25) 


C. Case II: The flat Gautreau-Tolman-Bondi background 

Gautreau considered the metric outside a massive body embedded in a universe containing inhomogeneous, spheri¬ 
cally symmetric dust and cosmological constant, A. Unlike the McVittie solution, the dust is in this case pressureless 
everywhere. The Gautreau solution is the flat limit of the Tolman-Bondi model written in curvature coordinates. It 
is also the nP = — A limit of the Leibovitz [46], solutions. Although not given here, we expect a treatment of dilaton 
evolution in asymptotically FRW, locally Schwarzschild, Leibovitz backgrounds to proceed along similar lines to this 
case. In co-moving coordinates, which are the most appropriate for considering the exterior expansion, the metric can 
be written as 


ds 2 = At 2 - R 2 r {t, r)dr 2 - R 2 (t , r){d<9 2 + sin 2 ddcj) 2 }, 

where 


R\ = 


2m + 2Z (r) 
R 


^AI? 2 . 

3 


The energy density of the dust is given by: 


ne = 2 Z tr /(R 2 R tr ). 


(26) 


(27) 


(28) 


We use the remaining freedom we have in the definition of R to prescribe that at some epoch of interest, t = to, r 
is the physical radial coordinate, so R(to,r) = r and the surface (t = to,r = const ) has area 47rr 2 . The surface of 
our spherical massive body is at R = R s . The quantity m + Z (r) is then defined to be the active gravitational mass 
contained inside the shell (t = to, r) = const, and in is the active gravitational mass of the central massive body at 
t = to- In general, the mass of the central object will change over time due to accretion of external material. We 
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require that as t —* oo the central mass remains strictly positive. The dust density must also be everywhere positive 
and tend to homogeneity for large r. Therefore we have the conditions: 


Z{R S ) 

= 0 

(29) 

Z r 

> 0 

(30) 

lim Z(r) 

r—* oo 

- if> (0) H 2 r 3 

~ o ^ “dust^O ' 

(31) 

lim Z(r) 

> — m 

(32) 


The last requirement is only strictly necessary if we wish the interior to have the required properties for t^> to- All we 
actually require for this analysis to hold is that the interior remains Schwarzschild to leading order over the time-scales 
appropriate to our intermediate matching region. Under our assumptions, equation (27) has exact solution: 


R(t,r) 
ti{r) 


6 (m + Z) . 2 ( v / 3A 


1/3 


to — 


A 

2 

V3A 


sinh 


-(t — h(r)) 


sinh 


r 3 A 


6 (m + Z) 


(33) 

(34) 


1. Exterior Metric 

As in the McVittie case, we define dimensionless coordinates p = H^r and r = Hot,. With respect to to these 
coordinates we can then asymptotically expand Z(r), order by order, in the small parameter <5. First, we write: 

H 0 Z(r) ~ \^ st p 3 + S p Zi(p) + o (S p ). 

Next, we require 2S p zi(p)/Q < '^ l l st p 3 ^ 1 for P ~ 0(1) so that this is valid asymptotic expansion; this ensures p > 0. 
The unperturbed spacetime is then FRW. Within the framework of a given model it will be possible to find the value 
of p. Given the expansion of Z we can then expand R order-by-order using equation (33). Putting this expansion 
back into the metric, eq. (26), we will have the expanded exterior metric in the form: 

d 4xt ~ Ho 2 [g^b {r, p) + 5 p g { ab (t, p) + o{5 p )^j dx a dx b , 


where g^ 


is the FRW metric. 


2. Interior Metric 

We define T = (t — t 0 )/Lj, and £ = R/R s to be our coordinates in the interior. We express Z(r) = Z(£,T) and 
expand out in power of the small parameter <5: 

Z(r)/m ~ c>Vi(x) + 0 (<5 9 ) • 

Where, from eq. (33), x = (£ 3 ^ 2 — 3T/2) 2 ^ 3 ; R s \ = r + 0(S q ). Putting this expansion back into the metric, eq. (26) 
gives the first two terms in the asymptotic expansion of the interior metric: 

dsL ~ R l [jab (£) + 6q Jab (£» X) + o(^ 9 )) dx a dx b , (35) 

where jty is the Schwarzschild metric with mass in/R s in Painleve-Gullstrand coordinates (with a re-scaling of the 
time coordinate): 


j i a b ) dx a dx b = ^-d T 2 - (d£ - C 1/2 dT^ - $ 2 {d(9 2 + sin 2 6»d0 2 }, 
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and Jab is g iven b y : 

dx a dx b = ^dCdT-^MdT 2 . 

Unlike in the McVittie solution, the next-to-leading order term in the interior metric here is not quasi-static, however 
by defining T = sjR s /2mT to be the time coordinate, the metric can be seen to be quasi-static in a different sense. 
Instead of only depending on T through ST, it depends on T only through the combination 2 mT/R s , and for most 
bodies, with the exception of black holes and neutron stars, 2 m/R s <C 1. For eq. (35) to be a valid asymptotic 
expansion we require S q \i i(%) <C 1 for % ~ 0(1). 


3. Energy - Momentum Tensor 


The Gautreau model is a subcase of the Tolman-Bondi dust models and so the pressure vanishes, P = 0, and 
= e. In the interior it is given by: 


R 2 kT^ = S q e\^, X )+o(S q ) (36) 

{2mS q \ Mx) lX 

= \~ii7) 

An expression for the exterior expansion of is not given here, but it can be found, with reference to the exact 
solutions for R(r,t), in the same way as the exterior expansion for the metric. 


D. Case III: Tolman-Bondi models with simultaneous big bang 


The final class of specific model that we will consider in this paper is, like the Gautreau solution, also a subcase 
of the Tolman-Bondi models [43]. This class is distinct from the Gautreau models in that we demand the big-bang 
to be simultaneous for all observers, and we do not require the hyper-surfaces of constant t to be flat. In the limit 
of spatial homogeneity these models can therefore reproduce all dust-plus-A FRW models. In the Gautreau model 
-R t > 0 everywhere. This implies that the world lines of the matter particles will stream out of the surface R = R s \ 
swept along by the cosmic expansion. In contrast, in this model we can, and will, require R t < 0 near R = R s . As a 
result, dust particles will fall onto our massive body; this seems a more physical reasonable scenario than that of the 
Gautreau case. In co-moving coordinates the metric can be written as 

ds 2 = dt 2 — —^-dr 2 — R 2 (t, r){d# 2 + sin 2 0d</> 2 }, 

1 — k{r) 

where 


s 2 , = -Mr) + — 


2Z(r) 


R 


5 “ 2 - 


The matter content of these models is pressureless dust with a cosmological constant. The dust energy density is, 
as in the Gautreau models, given by equation (27). As in the previous example we use the freedom we have in the 
definition of R to prescribe that at some epoch of interest, t = to, R is the physical radial coordinate i.e. R(to, r) = r 
and the surface (t, R) = const always has area 47 tR 2 ; this requirement, combined with the simultaneity of the initial 
big-bang curvature singularity , determines the form of k(r). The surface of our massive spherical body is at R = R s . 
As before, to + Z(r) is the active gravitational mass interior to the surface (t = to, r) = const and to is the active 
gravitational mass of the central body at f = to- The mass of the central object will grow over time as a result of 
accretion. We also require the dust density to be everywhere positive and tend to spatial homogeneity for large r, 
hence we need 


Z(R S ) 
Z r 

lim Z(r) 




0 , 

0 , 



(o) 

dust 


H 2 r 3 . 


(37) 

(38) 

(39) 


We want the zerotlr-order, exterior limit to be a FRW spacetime with curvature parameter k; this requires 
linir—xx, k(r) = kr 2 . We must require k(r) > 0 in the interior region; however we do not require k(r ) > 0 every¬ 
where. The exact solution for R(r,t) was found by Barrow and Stein-Schabes and is given in ref. [47]. 
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1. Exterior Metric 

As in the previous examples, we define dimensionless coordinates p = Hqt and r = Hot, and expand Z(r) order by 
order in the small parameter 5. We write: 

H 0 Z(r)^lil^ + SP Zl (p)+o(^). 

We require 26 p zi(p)/fl^ st p 3 <C 1 for p ~ 0(1) to ensure that this is valid asymptotic expansion; i.e. p > 0. The 
unperturbed spacetime is then FRW. The value of p is model dependent. With Z(r) specified, we can use the exact 
solutions of Barrow and Stein-Schabes, [47, 48], for R(t, r), and the requirement that R(r, to) = r to find the expansion 
of k(r) and from there the expansion of the metric; schematically we have: 

k(r) ~ n^p 2 + 6 p E 1 (p)+o(6 p ), (40) 

ds 2 ext ~ Hq 2 (g^ (r, p) + 8 p g [^ (r, p) + o(S p )^j dx a dx\ (41) 

where is the FRW metric. 

2. Interior Metric 

We take T = (t — t-o)/Li, and £ = R/R s to be our coordinates in the interior, and: 

Z(r)/m ~ S q pi(g) + o(S q ), 

where g = (£ 3 / 2 + 3T/2)“^ 3 ; R s g = r + 0(8 q , 6 2 / 3 ). From the exact solutions we find: 

k(g) = 5 2 / 3 k 0 (1 + 6 q p 1 (g) + o ( S q )) + O (<S 5/3 ) , 

where 



As with the Gautreau case, j^ is the Schwarzschild metric in Painleve-Gullstrand coordinates. The discussion about 

the quasi-static nature of the Gautreau term, given at the end of section IV C2, applies equally well here. For 
the approximation above to be a valid asymptotic expansion we require 5 q p\(x) ^ 1 for X ~ 0(1). 
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3. Energy - Momentum Tensor 

The energy density of the dust locally is given by: 

R 2 s*t£ = 


We will not give the exterior expansion of Tff here explicitly, although it can in principle be found with reference to 
the exact solutions for R(r , t) if required. 


S q 4&x)+o(6«) 
( 2 mS q \ AT ix), x 
) e / 2 x i/2 


(44) 


o(6«). 


E. Boundary Conditions 


As the physical radius tends to infinity, R —> oo, we demand that the dilaton tends to its homogeneous cosmological 
value: (f>(R,t) —> 0 c (f). We can apply this boundary condition in the exterior region but not in the interior. To solve 
the interior dilaton field equations we need to specify the dilaton-flux passing out from the surface of our ‘star’ at 
R = R s . At leading order we parametrise this by: 

-R 2 S ^1 - d fl 0 o | R=Rs = 2 mF (0 O ) = ^ ' dR'R' 2 B t< t,(MR'))^(R'), (45) 

where 0o = 4>o(R = R s )- The function A(0) can be found by solving the dilaton field equations to leading order 
in the R < R s region. If the interior region is a black-hole ( R s = 2in) then we must have F((j>) = 0; otherwise we 
expect F((j>) ~ I? ^(0). Without considering the sub-leading order dilaton evolution inside our ‘star’, i.e. at R < R s , 
we cannot rigorously specify any boundary conditions beyond leading order. Despite this, we can guess at a general 
boundary condition by perturbing eq. (45): 


-R 


2 

s 



9rS(^) 




5 (V—gg RR ) 9 r4 >o 


+ 2 S(M)F (<^>o) + 2toI 7, j 0((^o)^ (</>o) + smaller terms, (46) 


where S(X) is the first sub-leading order term in the interior expansion of A; M is the total mass contained inside 
£ < R s and is found by requiring the conservation of energy; at t = to we have M = m. Only 5 (cf) o) remains unknown, 
however we shall assume it to be the same order as 6(<j)) and see that this unknown term is usually suppressed by a 
factor of 2 m/R s relative to the other terms in eq. (46). 


V. APPLICATIONS 


A. Zeroth-order solutions 

In accord with our prescription, all of the models that we have considered share the property that, to lowest order 
in S, the interior is pure Schwarzschild, and the exterior is pure FRW. 

1. Exterior 

In the exterior, the dilaton field takes its cosmological value: (f> = <f>(t) c to zeroth order and so satisfies the 
homogeneous conservation equation: 

d 2 (f)(t) c + 3Hdt(j)(t)c = B ( 0 C ) Ke c dust (t) - V $ ( 0 C ) , 

The effect of the interior region on the exterior should, even for finite 5, become increasingly smaller as r —-> oo. As 
a result there will be no sub-leading order, r-independent, terms in the exterior expansion of 0. Equivalently, the 
homogeneous mode of 0 in the exterior will be given by the cosmological term, 0 C , alone to all orders 
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2. Interior 

In the interior the dilaton field satisfies the wave equation on a Schwarzschild background at zeroth order in 5; we 
take <j> to be quasi-static to leading order. By applying the zeroth-order boundary condition, eq. (45), we find: 

4>o = </>e (ST) - F (<f> o) In ^1 - , 

where (j) e (ST) is a constant of integration to be determined by the matching procedure. 

3. Matching 

As specified above, we should rewrite the exterior and interior expansion of </> hr some intermediate region, with 
length scale, Lj < L int < Le ■ We define 

LintT = (t — to). 

The only homogeneous mode in the exterior expansion was the cosmological term, (j> c (t ) = 4>hom (LintT/L e). This 
term will therefore appear at leading order in the intermediate expansion. Any other homogeneous terms that result 
from taking the intermediate limit of the exterior expansion must be sub-leading order. Therefore, we conclude that 
4>hom = 4> c (t) is the only leading order, homogeneous term in the intermediate expansion of 4>. In addition, all sub¬ 
leading-order homogeneous terms must also depend on time only t , and so will be quasi-static (Li nt T/L e dependent) 
in the intermediate regions. 

When the intermediate limit of the interior approximation is taken, it is clear that homogeneous terms in the 
interior will map to homogeneous terms in the intermediate zone. The matching criteria therefore implies that all 
homogeneous terms in the interior must be quasi-static, i.e. only depending on time through ST, and that at leading 
order: 

4> e (ST) = ct> e (L int T/L e ) = 4>hom (L int T/L E ) = 4> c (t) 

The interior solution therefore reads: 

1 So = <t>c (t) - F (4>o) In ^1 - , 

and 4>o = 4>c ( t) — F (<fo) In ^1 — . We see directly the effect of the cosmological evolution of <j) on the local region. 


B. Case I: The McVittie Background 

When the matched asymptotic expansion method is applied to the McVittie background, the analysis goes through 
in much the same way as it did in example 2 above. In the interior we find: 

^ = 4 o ) +<^#+ 0 (< 5 4 ), 

lbs 

where (ft'p satisfies: 


A 


(i)\ _ />e (ST) 


2m h’(ST) 

'ic s AV*m 


h (ST) 


3 | 1 2 m 

A 1 / 2 (0 + 2f? s AA 3 /2(£) 


4>e (ST) 




F(4> o) 


, 4 , (4° } ) 12ft2 + 6 h’Ad) 1 ^ + 3 n k (ST) 


A(OJ 

1 + 3AV 2 \ (IF A V 2 \ 5 


- 12Q a 


+R 2 s V, 4 , ($) . 
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This can be integrated once the functions B ^ and V t $ are specified. If we also have 

B^E)F(t o)f| 

b^(4>e) ^ 

- —— <sr 1 

vA4>e) ’ 


(48) 

(49) 


then we can solve for <j>j as an asymptotic series in 2 m/R s (provided this quantity is small). These relations will 
almost always hold provided that 2 m/R s £, <C 1 and the cosmological value of </> does not lie near the minimum of V 
or B. To O (2m/R s ) we find: 



+\bAM ('•'(«’) + x 



( C(6T) 


D(ST) 


h'F (0o)) + \V k (8T)F (0o) - /i 2 .F (4) 

^k(8T)) — —F (fo) B^ (4 >e) F s K ^du St + F (^o) R s V,<t><t> 



(50) 


where C(ST) is determined by the boundary conditions on R = R s and D(8T) comes from the matching conditions. 
The matching procedure would set D(8T) = 0 if we were to continue our asymptotic expansions to an order smaller 
than we do here. We do not expect C(ST ) to be any larger than the other terms in this expression when £ = 1. The 
exterior expansion, to 0(8), along with the matching conditions that determine the otherwise unknown functions in it, 
is given in Appendix A. In this appendix, we show that these matching conditions are self-consistent, and hence that, 
to the order considered, the matching procedure works correctly. As with the simple example we considered in section 
HID 2, we see that when the matching conditions are applied, the background sets the form of the homogeneous terms 
in the interior solution, and the interior solution gives us the behaviour of the inhomogeneous terms in the exterior 
expansion. This is to be expected, since in these models the interior region is the sole source of the inhomogeneities 
in both the spacetime and the dilaton field. This picture would be more complicated if we were to consider more than 
one interior region, as in the case of a nested series of shells of matter of differing densities. 

Our original concern was the behaviour of the time derivative of </> in the interior and we have shown that in the 
McVittie background: 

~ 1 — F >4> (0o) ln(l - 2m/R) + O (s 2 (2m/I?)) 2 ) . 

Since, in most cases of interest, the dilaton to matter coupling is weak and we are far from the Schwarzschild radius 
of our ‘star’, both F((/j> o) and 2m/ R are <C 1. Hence: dt<f>i / dt<f> c ~ 1, and the local time evolution of the dilaton 
tracks the cosmological one. The strength of this result arises partly from our restrictive choice of background metric. 
If the background is taken to be Tolman-Bondi rather than McVittie, we can find quite different behaviour. 


C. Case II: The Gautreau-Tolman-Bondi Background 


We assume that we have specified some Z(r; 8) that has outer and inner approximations as given in sections IV C 1 
and IV C 2 respectively. But we must be careful to ensure that the form of Z(r; 8) is such that these two approximations 
can be matched in some intermediate region. We can now proceed to solve the dilaton evolution equation in the inner 
and outer limits and then apply our matching procedure. We have already seen that at zeroth order the matching 
conditions imply that <f> is approximated by 


inner approx ~ (£; <5T) + O (5 q ) 

outer approx ~ (j) c (r) + O (8 P ). 


<t>c (t 


t 0 + ST) - F (<f>o) In 


2 m 

R4 


+ 0 ( 6 *), 


We now consider the next-to-leading order terms with a view to determining how time variation of (f> in the interior 
is related to its cosmological rate of change. 

In the interior the local energy density and metric are, to O (S q ), depend on T only through \ and are written as 
functions of £ and The calculation is easiest if we move from (T, £) coordinates to (T,\) ones and write: 

0 ~ 4 0) (6 ST) + 8^P (T, X ; ST) + o (8 *), 
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where <j>] satisfies (for q < 2): 


2 TO 

i?7 



1 / 

'« 6/ V> s 

k 

II 

to 

%. 

'eT 

Mi (x), x 

x l/2 \ 

.x 1 / 2 ^ 

x 1 / 2 

( 2 m N 

2 

) F(j 0 ) 

Mi (x), x 

2mi (x) 

\ F* / 

(i - tt) 

« s/2 ( 1 --©) 2 


(51) 

(52) 


There is also a single term that appears at order S which we have omitted from the above expression; we shall deal 
with this later. In general, the above equation is difficult to solve exactly, however, in most of the cases of interest 
the surface of our ‘star’ will be far outside its Schwarzschild radius and so 2 m/R s <C 1. We may then solve eq. (51) 
as an asymptotic series in 2 m/R s , valid where £ Given the nature of our problem we are only interested in 

the leading term: 



2 to „ . , . 

(4>c) 

rC s 



Mi (x') lX 
ax\T) 


Mi (x) 
£ 


+ C(T)j + O 



Here, A(T) and B(T) are ‘constants’ of integration; A(T) should be determined by a boundary condition on the 
surface at R = R s i.e. at £ = 1, and B(T) will come out of the matching. We argued in section V A that the matching 
requires that all homogeneous terms in the interior expansion be quasi-static, thus B(T) = B(ST). Without knowing 
more about the region £ < 1, we do not have a boundary condition capable of determining A(T). If, however, we 
assume that the prescription given by eq. (46) is at least approximately correct then we can proceed. We assume that 
the next-to-leading order perturbation of <j >o occurs at the same order in 8 and 2 m/R s as the next-to-leading order 
term in <f>. Given these assumptions we prescribe our full boundary condition at R = R s to be given by eq. (45) but 
with (f>o —> jo + O (S q 2m/R s ) and to —> m + Z\r=r b . The resultant boundary condition on (jyp is then remarkably 
simple: 

(x(£ = 1 ,T)). 

With this choice of boundary condition, 


A(T) = (1 - F(j 0 )/B(<f> c )) Ml (x(£ = 1 ,T) 


and the interior solution is: 



-$^,*(0 C ) 

ILa 



d A 

* Z(X',T) 


Mi (x) 


F(jo) \ Ml (x(g = l ,T) + h{8T) 


+ o 



The unknown quasi-static term, B (ST), will not affect the leading-order time-variation of <j>P and so will not have 
any bearing on our result. We could, in principle, solve for </> in the outer approximation, however, since we are 
interested only in dynamics of the dilaton held in the interior, but we will not do so. 

Recall that we dropped an order-5 term in writing down eq. (51); this perturbation decouples from the one at order 
S p , and is quasi-static. Writing 


we find: 


<h ~ 4>P(ST) + + SpPit, ST) 

P fa - 2 m/R s )pp) ? (Ve m - (fo) In 




Solving this we find: 
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f 2m\ 
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j'e (ST) 
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(53) 
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where A is a constant of integration and 


m 





/ —2 



2 to 


3/2 


In 



2 


2r?z 


1/2 

— 2 In 



The term proportional to l(£f) dies off as , and is suppressed by a factor of 2 m/R s ^ relative to the first term in 
eq. (53). Both these terms only depend on time through ST and are thus deemed quasi-static. These terms will not , 
therefore, change the leading-order behaviour of the time derivative of <j> in the interior. 

Although we are concerned mostly with objects which are much larger than their Schwarzschild radii it would be 
nice to be able to address the problem of black-hole gravitational memory via this method. For the zeroth-order 
approximation to <f> to be well-defined on the horizon we need F(<j> o) = 0 (this is just a statement of the ‘no-hair’ 
theorem for Schwarzschild black holes). We can remove any divergence in near the horizon by an appropriate 
choice of the constant of integration, A. We take A = cf>' e (ST). This is analogous to what was done by Jacobson in 
[35]. We can now absorb the first and last terms on the RHS of eq. 53 into the definition of f> e {5T) by a definition of 
the time coordinate: 


where 


4>e{ST) -► MST) 


T = T+ 2 





The zeroth-order matching is not affected since S(T — T) will remain sub-leading order in any intermediate region. 
We note that LjT = v — R—2m ln(i?/2m) where v = t s + R + 2mln(/?/2m — 1) and t s is the usual Schwarzschild time 
coordinate. Thus the leading-order homogeneous term in </>/ has the same behaviour as that predicted by Jacobson 
(see ref. [35] and section IIB above). 

By expanding <j) e w.r.t. the co-moving time, LjT , rather than the Schwarzschild time, t s , as Jacobson did, we avoid 
one of the problems encountered in his analysis: that the time coordinate use to make the expansion diverges on the 
horizon. 

At lower orders we can always ensure, by choice of the constants of integration, that lim/j^ 0 (l ~ 2m/R)dR<j>\T = 0. 
This boundary condition ensures that the inner approximation does not diverge as R —+ 2to, and remains valid on the 
horizon. This justifies the statement that, if p > 1, Jacobson’s prediction is itself a valid asymptotic approximation to 
4> on the horizon. The condition p > 1 is equivalent to: e; oca ;/e c <C l/(2mHo)(f$> l)i with e; being the average value 
of the energy density, outside the black-hole, over length scales 0{2m) from the horizon. 

We are now in a position to study the conditions under which the local time variation of <f> tracks its cosmological 
value. In the interior we have: 


{f>i f 


W e (6T) (l-Fj (4> 0 ) hi (l ^|)) 


+S g B. 


2m 

’/A 


d.V 


e /2 (x r ,T) 


+ i 


F(M \ Mi,x(£=l ,T) 

b{ 4> c ) J e /2 


(54) 


The excluded terms are then certainly smaller than the second term, however, since they might he larger than the 
first term, both numerically and in the limit S —+ 0 , the above expression is not a formal asymptotic approximation. 
When the first term in eq. (54) dominates, condition (5) holds, and the local time evolution of the scalar held is, to 
leading order, the same as the cosmological one. Assuming our background choice is suitable for the application of 
this method, condition (5) fails to hold if, and only if, the second term in this expression dominates over the first. 


D. Case III: Tolman-Bondi models with simultaneous big bang 

In the interior, the leading-order corrections to the metric occur at either order S q or J 2 / 3 . The order S 2 / 3 correction 
comes from the leading, quasi-static, and ^-independent term in the expansion of /c(r); we have seen that this can be 
transformed away by a redefinition of the local time coordinate, T —> T*, as done in section IVD 2. 

In the intermediate region T ~ \/\ — 5 2 / 3 koT* , the zeroth-order matching goes through in the same way as it did 
it the two previous cases. The zeroth-order interior solution is 

</>o = <t>e (sVT^SV^oT*) - F (4) In ^1 - ■ 
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The analysis of the O ( 6 q ) terms then follows through in much the same way as for the Gautreau case, but with 

X = ^ 3/2 - 3T/2^ 2/3 -► T] = (C 3/2 + 3T/2y /3 . 

One hnds that the O ( S q ) correction to the interior solution is given by 


4 1} - jr B ,4>{<l>c) 

rC s 


r>_ 

dry 


Z(V',T) ? 

There is also a quasi-static order 6 correction to <f>i, just as for the Gautreau case. We find 


b { P = — v 7 1 — <5 2/3 fc 0 2 
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4>e (ST*) 
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3/2 
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R4 


with l(£) as given in section V C. As above, we shall choose the constant of integration, A , so as to make <f>f well-defined 
as R s l ; —> 2m whenever F (0o) = 0; we can then absorb this correction into the definition of <j> e : 

<f>e (Sy/l - S 2 / 3 k 0 T*^j - (f> e (Vl - A 2 / 3 fc 0 T*) . 


The transform T* —» T* is the same as the one for T —> T given in section VC; LjT = v — R — 2mln(i?/2m). As 
in previous case, we can now analyse the conditions under which the local time variation of (f> tracks its cosmological 
value: 


0/,f. « Sy/1 - SV 3 k 0 ct>' e (Vl - SV^koST*) ^1 - F t4> (0 O ) In (l - 




. 2m 

:) i?7 


/ F(0 o)A = 

7 C 5 / 2 ( ry ', T ) + V S ( 0 C ); ^/2 


+ ... 


(55) 


This not a formal asymptotic approximation, but the excluded terms are certainly smaller than at least one of the two 
terms in the above expansion. The effect of the transform T —> T* is to induce a slight, sub-leading order, lag in the 
time evolution of the scalar field, that increases as £ decreases. As in the Gautreau case above, condition (5) will fail 
to hold if, and only if, the second term in this expansion dominates. We will interpret this condition, and the similar 
one which arose from the analysis of the Gautreau case, in terms of what it requires for the local energy-density in 
section VII below. 


VI. CONDITIONS FOR THE ASYMPTOTIC EXPANSIONS 

By considering the growing modes in the interior expansion for <j), and the singular modes in the exterior expansion, 
we can say something about the position of the matching region, and in some cases show that a matching region could 
not exist. 


A. Case I: McVittie background 

At 0(5 2 ), the interior expansion of <fi in the McVittie metric has a mode that grows like £. This asymptotic 
expansion will certainly cease to be valid if, when rewritten in some intermediate scaling region, where £ oc 5~ a £, say 
with 0 < a < 1, terms appear that scale as inverse powers of S. In such a region 5 2 t; ~ <I 2-a £, and so such terms 
will always be suppressed by at least a single power of S. The zeroth-order term, <pf\ in the interior will itself have 
ceased to be a valid asymptotic approximation to <j> when the 5 2 £ terms dominate over the l/£ term. This will occur 
only when £ oc (V 1 , i.e. in the exterior region. At order A 4 , the fastest growing mode will go like £ 3 , and this will 
also only conflict with the first two terms in the expansion when £ oc <5 -1 . We conclude that eq. (50) is a valid 
asymptotic approximation to cf> everywhere outside the exterior region. Similarly we note that the most singular term 
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in the first two terms of the exterior expansion, as given in Appendix A, goes like 6/p. At O (d 2 ) we would find terms 
that behave as S 2 /p 2 . The exterior expansion will therefore only cease to be valid when p ~ 6, i.e. in the interior 
region. Therefore, in the McVittie case, the position of the intermediate region is not important; we can choose any 
intermediate scaling and the matching procedure will be valid. As a final check on our method, we have explicitly 
performed the matching for the McVittie background in Appendix A, and shown it to be self-consistent. 


B. Cases II and III: Tolman-Bondi models 

We assume that as pi (x) ~ \ n as X ~> 00 for some n > 0. At order 5 q , the growing mode in the interior 
approximation will grow like 5 q x n /£ in the Gautreau case, or as S q r] n /£ for the simultaneous big-bang models. 
Assuming that in some intermediate region 77 , x> £ ~ S~ a with 0 < a < 1 we can see that this growing mode will 
dominate over the zeroth order l/£ when a = q/n, and the asymptotic approximate will fail altogether if n > 1 and 
a > q/(n — 1). We did not explicitly find the exterior expansion of <fi, since we were only really concerned with its 
behaviour in the interior, however the first non-homogeneous mode should behave as 5 p z\(p)/p if p < 1 or S/p if 
p > 1. If zi(p) ~ p m as p —> 0, then the exterior expansion will break down if 

m < 1 — p a < 1 — —^—, (56) 

1 — m 

m > 1 — p a = 0. 

In addition to this consideration, we note that the transformation T —> T* used in case III will not be well-defined 
if £ oc S~ a where a > 2/3. Just by considering the behaviour of the next-to-leading order terms we can say that if 
both the interior and exterior zeroth-order approximations are to be simultaneously valid in some intermediate scaling 
region we need, for the Gautreau case: 


max 



< a < 


and for the simultaneous big-bang case: 

max ^0,1 — -j— ^^ < a < min ^—,2/3^ . 

If such an intermediate region does indeed exist then the zeroth-order matching performed in section V A will be 
valid. The general form of the interior approximation to O ( S p ) will then be correct; the only unknown function in 
this term is B{8T). If the matching works to order 5 P , as well as zeroth-order, then we have argued that B(ST) will 
be quasi-static. If the matching procedure does not work to this order then its quasi-static character may be lost. 
However, we would not expect it to vary in time any faster than the other O ( S p ) terms. So long as we can match 
the zeroth-order approximations in some region, we can find the circumstances under which condition (5) holds by 
comparing the sizes of the two terms in (54) and (55), for the Gautreau and simultaneous big bang cases respectively. 


VII. INTERPRETATION AND GENERALISATION 


We assume that 2mF(<f>o)/R s <C 1 in the Tolman-Bondi cases II and III considered above and that the matched 
asymptotic expansion method is valid. Condition (5) will then only fail to hold if the O (S q ) term in expressions (54) 
and (55) is larger that the order 5(j>^' term. For condition (5) to hold we require: 
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(57) 


where R S X , RsP ~ r and i? s £ = R\ t — to = LiT. For any given model this can be evaluated. As the expression is 
currently written its physical meaning is somewhat obscured. Firstly, the lower limit on the integral is not specified; 
we merely specified that the ‘constant’ of integration that arises from this lower limit should be quasi-static. We 
should therefore take the lower limit to be some value of r such that the contribution to the integral from that lower 
limit vanishes. If p > 1, the above criterion (57) will certainly hold. Therefore, in looking for where (5) fails to hold, 
we are concerned only with the cases where p < 1. In these situations, and indeed also for if 1 < 6 < 2, the local 
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matter density close to R = R s will be much larger than the cosmological dust density. We note that in the interior 
the energy density, e(r,t) and also energy overdensity, Ae (as p < 2), is given by 

K e(r, t) ~ = n(e - e c ) ~ 2 md p r3 ^ 1/2 - 

We also note that i?,* =~ ±(2m/R) 1 ^ 2 and R. r ~ r 1 / 2 /^ 1 / 2 . Therefore, in the interior region, we can substitute 
±(2m/i?.) 1 / 2 /xi, r for A(i? t «e) := — HRe c ). Doing this we find that our criterion, eq. (57), becomes: 

B A^c) (f dr'R^A( K R, t e(r',t)) + (l - Trjf$j) dtKAe) R=R 

-V- ; -V- •AAAI -f- L < i. (58) 

4 >c{t) 

The two expressions, (57) and (58), are equivalent whenever the interior approximation holds. These should be 
checked in the context of a given model. We can further simplify expression (58) by noting that, to zeroth order in 
2 m/R s , we should expect F(<j> o) ~ B^(0 C ), and so the second term in the denominator will be a factor of 2m/£s 
smaller than the first and so can safely be neglected. 

Expression (58) is especially useful since it is written in terms of the matter overdensity, ne. Using this expression 
wc will now conjecture a more general criterion, that will coincide with (57) and (58) whenever the matching procedure 
works, but which we believe should also apply to many cases where the matching would not have worked. Eq. (58) 
only features terms that are well defined for all Tolman-Bondi models, i.e. I? j7 , f,t and kAe. The condition 


I dr' R r A(R t ne(r', t)) < oo 

J OO 


(59) 


will hold provided kA e,t ~ o(l/i? 2 ) as R — » oo; this is certainly true whenever the extra mass contained in the 
over-dense interior region is finite. The integrals in (57) and (58) are specified by the condition that any constant 
term that remains as R —» oo should be neglected. So long as oc x" with n < 5/2, these integrals will be 
equivalent to the expression given in eq. (59), at least in the interior region. The analysis of section VIB suggests 
that if n > 5/2 and q < 1 the size of any potential intermediate matching region will be strongly bounded; thus, for 
most cases where we expect the matching procedure to work, n < 5/2. 

In the Tolman-Bondi models the radial velocity of the dust particles, as measured by an observer at fixed R, is 
given by R t ■ Assuming that the second term in denominator of (57) and (58) is negligible compared to the first, 
we can rewrite our criterion for condition (5) in the more concise and revealing form by replacing R t with the dust 
particle velocity, v: 


§£,»(&) |A(me)| ^ 

4>c{t) 

where dl = dr'R r and the path of integration, j(R), runs from R to oo, and A (vne) = vke — HRe c . Since we expect 
that (f> will only feel the effects of events that happened in its causal past, we should take 7 (i?) to run from R to 
spatial infinity along a past, radially directed light-ray. For the Tolman-Bondi models we have considered, the LHS 
of (60) will coincide (at least in magnitude) with the LHSs of (57) and (58). 

Eq. (60) has the form that we might expect for condition (5) to hold inside a local, spherically-symmetric, inho¬ 
mogeneous region produced by an embedded Schwarzschilcl mass in an asymptotically FRW universe: the LHS of eq. 
(60) vanishes as we approach the cosmological region and the local time variation of the dilaton field is driven by the 
time variation of the energy density along its past light-cone. 

We shall therefore conjecture that eq. (60) is a general sufficient condition for eq. (5) to hold: applicable even 
to spherically-symmetric (or approximately spherically-symmetric), dust plus cosmological constant, backgrounds in 
which the matching procedure would itself fail. 

If the cosmic evolution of the dilaton field is dominated by its coupling to matter, so \B^((j) c )Ke c \ » \V^(<j> c )\, then 
(j> c (t) = O (B^{4> c )Ke c (t) / H{t)) and condition (5) therefore holds, at the epoch t = to, whenever 

[ diff(i 0 )^ 3 «l. (61) 

e c\f Oj 

If, alternatively, the cosmic evolution is potential-driven, so \V^(j) c )\ \B^((f> c )Ke c \, then the LHS of the above 

expression will be suppressed by an additional factor of \B^((j) c )Ke c /V^((j) c )\ <C 1 . 
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In section IIA we reviewed Wetterich’s conclusion of [31] that condition (5) holds true when the cosmic evolution 
of the dilaton is potential dominated. Whilst we have identified gaps in the analysis performed in that paper, we have 
been able to establish a similar result. For a given evolution of the background matter density, we have seen that 
condition (5) is more likely to hold (or will hold more strongly) when \B^((f> c )Ke c /V ( f,((f> c )\ <C 1. The domination by 
the potential term in the cosmic evolution of the dilaton has a homogenising effect on the time variation of cf>. In the 
final section we shall apply eq. (61) to some physically relevant scenarios and show that, on the basis of the above 
analysis, we should expect </(x, t) ss 4> c {t) in the solar system. 


VIII. OBSERVATIONS IN OUR SOLAR SYSTEM 


We would like to apply out results to solve our original problem: whether or not condition (5) holds when the local 
region is the Earth or the solar system. We will consider a star (and associated planetary system) inside a galaxy that 
is itself embedded in a large galactic cluster. The cluster is assumed to have virialised and be of size R c i us t- Close 
to the edge of the cluster we allow for some dust to be unvirialised and still undergoing collapse. Since we have only 
performed our calculations for spherically-symmetric backgrounds, we should, strictly speaking, also require spherical 
symmetry about our star; for a realistic model this seems contrived. In an upcoming paper, [49], we shall consider 
more fully the effect that deviations from pure spherical symmetry have on our results; for now we assume that these 
effects are small, and we can relax the spherical symmetry requirement somewhat without invalidating our analysis. 
We consider the different contributions to the LHS of condition (61) in this astronomical set-up: 

I := [ dZ H(to) J 7T - T" = Iciust + I gal + Istar- (62) 

Jj(R) £ c(t o) 


For illustration we evaluate I for Brans-Dicke theory. For other theories, I will still take very similar values. We 
assume that the density of non-virialised matter, just outside the virialised cluster, is no greater than the average 
density of the cluster and if we move away from the edge of the virialised region, the over-density drops off quickly, 
i.e. as R~ s , s > 1/2. In this case, the Tolman-Bondi result, eq. (55), applies and the magnitude of I is bounded by: 


Rlust iS *0 1/2) \/ZMclust Rciust 


Cclust 


1.2 x 10 


-7 


( ^clust ^ 

f Vciust \ 

^ Mpc ) 

Vkms ” 1 / 


1.5 x 10 _3 [3/(2s - 1)]0“ 1 (1 + z vir f < 1, 


where we have used 3 GM / R c iust = v^. lust . R c i U st is the scale of the cluster post-virialisation, and v c i us t is the average 
velocity of the virialised dust particles in it; to = 13.7 Gyrs is the age of the universe. In the final approximation 
we used representative values R c i U st = lOOMpc and v c i us t = 200km s -1 appropriate for a cluster like Coma. Taking 
a cosmological density parameter equal to fl m = 0.27, in accordance with WMAP, and h = 0.71, we expect that 
for a typical cluster which virialised at a redshift z V i r <C 1, we would have I c iust ~ 5.7 x 10 -3 . The term in [] is 
unity when s = 2, i.e. 2GM/r —» const ; such a matter distribution is characteristic of dark matter halos. Different 
choices of s > 1/2 can be seen to only change this estimate by an 0{ 1) factor. As s —> 1/2 the matched asymptotic 
expansion method, and hence this particular evaluation, breaks down. We believe that the generalised formula for I 
will, however, still give accurate results. We have assumed that the baryon-to-dark matter ratio inside the cluster is 
the same as in the universe on average. If this is not the case, and the dilaton couples with vastly different strengths 
to baryons and dark matter, then I c i us t could be made as large as 3 x 10~ 2 , for these values of R c i U st and v c i uat ; this 
scenario would require the cluster to be comprised completely of baryons, and the dilaton to have zero coupling to 
dark matter. 

We now consider the galactic contribution. This will come about as a result of the galaxy slowly accreting matter 
from the intergalactic medium (IGM). We shall assume that the properties of IGM are the same as the average 
properties of the cluster. In this case the density of the IGM will be approximately constant and as such eq. (55) 
will not be strictly applicable. However, if we assume that our conjecture of the previous section holds then we can 
evaluate I ga i■ The IGM has an average particle velocity Vjgm\ as a result of this, the gravitational influence of our 

Galaxy upon it will only be felt significantly out to a radius Rq 11 ^ = 2 GM ga i/vj GM . We therefore take R = R,Q al \ 
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rather than oo, as the upper limit of the integral in eq. (61). 


l gal 


— 1 \J 2GM gal Rq 0,1 ^ -Z£M. 

-2.4xl0- 15 O- 1 ^(l + ^ r ) 3 

— 4.4 X 10 _6 (1 + Z v ir ) 3 <C 1, 


kms 1 

Vclust 


where we have taken (jgm ~ £ c iust and Vigm ~ v c i ust . In the last line we have taken M ga i = 1O 12 M 0 and v c i ust = 
200km s^ 1 as above. 

Finally we consider the contribution that results from our star accreting matter from the interstellar medium (ISM). 
This calculation proceeds in the same way as the one for I ga i- 


—2t^ 1 \/2GmrG——— 

£c 

— 1.2 x lO” 13 !!” 1 -^- 
Mg 

-1.8 x 10“ 14 n< 1, 


kms 


-l 


vism 


h 2 n 


where tq = 2 Gm/vj SM , and €ism = n protonscm -3 where n ~ 1 — 10 4 . We have taken vjsm = 5kms _1 and 
m = Mq to give the final approximation. It is clear that, in general, I star <C I ga i Idust- The infall of dust into 
cluster will tend to be the dominant contribution to the LHS of eq. (61). The estimate of I c iust given above should 
be viewed as an upper bound on its value; even still we have seen that is small compared with 1, and hence that we 
should expect condition (5) to hold near our star; with deviations from that behaviour bounded by the 0.6% level 
(if the dilaton couples only to baryonic matter this could rise to as high as 3%). Assuming that the conditions in 
our solar system are not too different from those considered above, we conclude that irrespective of the value of the 
dilaton-to-matter coupling, and what dominates the cosmic dilaton evolution, that 


0(x,f) « </> c (i) 


will hold in the solar system in general, and on Earth in particular. We have seen that if the cosmic evolution of <p 
is potential dominated then this only serves to strengthen this result. The only caveat is that since we have ignored 
the back reaction of (f> on the background cosmological geometry we are limited to \B^\ <C 1. If B^, and hence the 
back reaction, is large then the scope for violations of condition (5) increases; if the coupling to gravity of the dilaton 
is large, then dilaton field may itself undergo gravitational collapse. 

We can also conclude, from the above analysis, that, even if condition (5) is violated during the collapse of an 
overdense region of matter, once the region stops collapsing as a result of its virialisation, the time evolution of the 
dilaton field tends towards homogeneity. Contrary to what has been claimed before in the literature, virialisation 
does not stabilise the value of the dilaton, and protect it from any cosmological variation. 

In addition whilst we did not explicitly calculate the interior solution for 2 m/R s = 0(1), the evolution eq. (51) is 
still valid in this case. We should therefore expect the magnitude of the rate of time variation in <f>j to be similar 
to what was found for the 2 m/R s <C 1 case. We therefore conclude that over time-scales that are large compared to 
2 Gm, condition (5) will hold whenever condition (61) does, and there will be no significant gravitational memory of 
the value of G from the epoch when black holes or gravitationally bound structures first formed. This result agrees 
with the prediction made by Jacobson, [35], the inhomogeneous models of [33], and the numerical calculations of 
Harada et. ah, [34] , who also studied the Tolman-Bondi background. Indeed we have seen that whenever the, not 
very restrictive, condition, e;/e c <C l/(2mffo)(> 1), holds then eq. (8), gives a true asymptotic approximation to the 
behaviour of the dilaton close to the horizon. 


IX. SUMMARY 

In this paper we have considered the extent to which a cosmological time variation in a scalar field (dilaton) would 
be detectable on the surface of gravitational-bound systems that are otherwise disconnected from changes that occur 
over cosmological scales. This problem is of particular relevance when the dilaton defines the local value of one of the 
traditional ‘constants’ of Nature, and when the system is the Earth or our solar system. Several scalar-tensor theories 
have already been developed which self-consistently describe the variations of traditional ‘constants’ of Nature, like 
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a and G. By matching rigorously constructed asymptotic expansions of the associated scalar field found in different 
limits, close to the Earth and far away from it at cosmological scales, we have been able to derive approximate, 
analytical expressions for the scalar field near the surface of such bound systems. This result was found under two 
major assumptions: the physically realistic condition that the scalar field should be weakly coupled to matter and 
gravity(in effect the variations of ‘constants’ on large scales occur more slowly than the universe is expanding ) and 
thus have a negligible back-reaction on the cosmological background, and the less realistic one of spherical symmetry. 
We do not expect the relaxation of the spherical symmetry condition to greatly alter the qualitative nature of our 
conclusions, and we shall present a rigorous treatment of the problem when no symmetry is present in a subsequent 
work. 

Finally we have extracted from our analysis a sufficient condition for the local time-variation of a scalar field, or 
varying physical ‘constant’, to track the cosmological one, and we have proposed a generalisation of this condition 
that is applicable to scenarios more general than those explicitly considered here. By evaluating the condition for an 
astronomical scenario similar to the one appropriate for our solar system, we have concluded that almost irrespective 
of the form of the dilaton-to-matter, and the form of the dilaton self-interaction, its time variation in the solar 
system will track the cosmological one. We have therefore provided a general proof of what was previously merely 
assumed: that terrestrial and solar system based observations can legitimately, be used to constrain the cosmological 
time variation of supposed ‘constants’ of Nature. 
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APPENDIX A: EXPLICIT MATCHING FOR THE MCVITTIE BACKGROUND 


1. Exterior Solution 

In the exterior region <f>E ~ <t >\'e + 54> { e + 0(S 2 ) where: 

^3 drfa^r^) - ~ d P ( , , l^n , d P^ 


a*p 

3/2 


l + lQ0 p 2 


= (I 7 ) ^ + &(r)p ^ — {^e" + H { e' + 3B 0 (^) (/i' + ^fifc(-r))} 


+B. 


E^Ktdust^E ~ V^i^E^E 1 - 


aOKA 1 ) 


Despite the complexity of this formula we can solve it admits separable solutions in r and p\ 

, (1) _/2mV' a (l + W) i ' 2 , 


RsJ 


-T(t) + $(rp). 


where T(r) satisfies the ODE: 

1 


d T (a 3 d T T) - ^^T-B^(4 o) )7?- 2 K el st T + H- 2 U 00 (4 o) )T 

U, 

= ^e' + He ] ' + 3 B^((j> E ) (ti + ^-Dfc(r)) | . 

When k > 0, $ = y^c 7 • T+(r)A+(p), c 7 G <C 2 : 


X 1 = 


pi'ya(p) 


sin a 


, x . -l( V^kP 

a(p) = sin 1 


1 + I^kP 2 


(Al) 


and T+ = (T+’ 1 , T+’ 2 /, where T+’ 1 and T+’ 2 are linearly independent solutions of the following ODE: 
±d T (a 3 d T T+’\T)) + ( 7 2 ^ 1 ^(rjT+’V) = 
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If k = 0. then <J) = djc(j) • T°(r)X°(p), c(y) G C 2 , and X°(p) = e llp /^p. T°(r) = (T° a , T°’ 2 )‘ satisfies: 

^d T (a 3 4 5 6 7 8 9 10 d T T® ,i ( t)) + 7 2 T 7 °’ i (r) = - H^VM^T^r), 

where i = 1, 2. Finally if fc < 0 then $ = dy \J — fl°c( 7 ) • T ~(t)X~( p), 0 ( 7 ) G C 2 : 

i'yrv 

X~{p) = - 


a = sinh 


sinh a 

1 (JEEp 


1+3 ^° k p 2 r 


and T-(r) = (T-> 1 ,T"> 2 ) t : 

(d 3 d T T~^r)) - ( 7 2 + 1) n k (r)T~ p (T) = B^^H^^T-^t) - 
where * = 1,2. In all cases we shall fix our definition of the T 7 by the normalisation: T 7 (tq) = 1. 




2. Matching Conditions 


By making our interior and exterior solutions for </> we see that: c 7 = c _ 7 for k > 0, and 0 ( 7 ) = c(— 7 ); z is 
the complex conjugate of z. Defining A n (r) = 7 n c 7 • T+ or A n (r) = f dj j n c(j) ■ T?/ for k > 0 and k < 0 
respectively we see that: 


Mr) = EE}. _ T(t) 
a 

Ai{t) = 0 <^=> c(y) = c(—7) 

fifc| ^2(r) = O' + 2 h 2 )aF$ 0 ) + ~ K00) «+(^o) + n k (r)aF(<j> 0) 

-5 (4 0) " + *</4 0) ' + 3^.0 (/i' + ^fifc(r))) - ifigT. 


(A2) 

(A3) 

(A4) 


In principle we can invert the equation for Aq(t) to find the c)( 7 ). 

Given the definition of the A n , the expressions for Aq and A 2 must satisfy a consistency relation (if there were not 
to satisfy this, the matching procedure would be invalid). We have checked that this relation does indeed hold here. 
The relation is: 


l^fcO)! A 2 - flk(r)A 0 = (a 3 d T A 0 ) + (BM^E^^dust ~ v ,</><!>) A o- 

We can see, explicitly, that the matched asymptotic expansion method works for the McVittie background. 
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